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ABSTRACT 


The head-flow rate curve of the partial emission pump is flat and with a sudden drop-off near the maximum flow rate. Thus, it 
is difficult to fit the experimental head-flow rate curve of the pump by using a polynomial, which is suitable for centrifugal pumps. 
To tackle this problem, a new mathematical model was put forward in the inspiration of the mathematical model of the equivalent 
electric circuit of photovoltaic cells/modules. The physics behind the model was described by using flow rate saturation effect in 
partial emission pumps. The model was implemented in Excel and applied to the experimental head-flow rate data of seven partial 
emission pumps with different specific speeds in a range of 17-80. The model can produce high-quality curve fitting (coefficient 
of determination of greater than 0.9), The model is adaptive to the variable experimental head-flow rate data of partial emission 
pumps due to changes in impeller structure, rotational speed, number of blades and impeller diameter. Further, the model can fit 
the experimental head-flow rate data with drooping effect at low flow rate in partial emission pumps satisfactorily (coefficient of 
determination of greater than 0.94). The model potentially can be used in experimental data processing in the pump industry and 
pump-system modelling for partial emission pumps in future. 
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1. Introduction dramatic sharp cut-off at the maximum flow rate in the experiment, 
as shown in Figure 2. Usually, a 2"*- or 3"-order polynomial, which 
is mostly applicable to best fit the head-flow rate curve of centrifugal 
pump, is invalid for a partial emission pump. How to fit this kind of 
head-flow rate curve is concerned in the article. 


Partial emission pumps, as shown in Figure 1, are special centrif- 
ugal pumps developed in the 1960s, and have found extensive appli- 
cations in aerospace, chemical, petroleum-chemical industries [1, 2], 
offshore technology [3] and so on. The pumps have an open impeller 
with straight radial blades in 90° blade discharge angle (measured 700 ; 
from the reverse direction of the impeller rotation) [4]. j 
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Figure 2. Experimental head-flow rate curves at the rotational speeds of 


Figure 1. Impeller and volute configuration of a partial emission pump, (a) : 
15000, 20000, 25000 and 30000r/min, the data are after [1] 


axial-cross-sectional view, (b) cross-sectional view in volute mid-span 


As a result, the pumps often have a flat head-flow rate curve with 
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2. Model and method 
2.1 Current-voltage model for photovoltaic cells or modules 


The author had conducted a few modelling studies on the electric 
circuit of flat photovoltaic (PV) cells or modules [5-8]. A flat mono- 
crystalline PV cell or module can be presented by using a single di- 
ode equivalent circuit with five combined parameters such as photo- 
current, I,,, diode reversal saturation current, Ig, diode quality fac- 
tor, M, combined series resistance, R,, and shunt resistance, Rs, 
[5], as shown in Figure 3a. The performance curve, i.e. [-V curve 
for the PV cell can be expressed mathematically by [5]: 

1 = Ipn — la [exp {2} — 1] - SR (1) 
Where V and / are the output voltage and current of the PV 
cell/module, respectively; T is cell temperature, g is the electron 
charge, g=1.60217646x10"°C, and k isthe Boltzmann constant, k 
=1.38065031x10J/K. Ip, depends on both solar radiation inten- 
sity and cell temperature in the silicon layer, while Ig is only cell- 
temperature dependent. A set of I-V dataofaPV module was fitted 
by using Eq. (1), and good agreement between the data and the 
model is achieved, as shown in Figure 3b. This fact suggests that the 
exponential term in Eq. (1) can properly describe the rapid drop in 
the current when the voltage is beyond a certain value. 
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Figure 3. Single diode equivalent circuit of a flat monocrystalline PV module 
(a), and the experimental current-voltage ([-V) curve of the module at the 
standard solar radiation intensity S=1000W/m? (b), the plots are adopted 
from [5] 


Physically, active electrons can be activated in a PV cell of semi- 
conductor materials at a solar radiation intensity and cell temperature. 
The active electrons can travel to the p-side from the n-side in the 
cell to form a current inside the cell. If the cell is connected to a cir- 
cuit, the active electrons will go back to the n-side to combine with 
the holes left by the electrons through the circuit outside the cell, thus 


a current can be established in the circuit. The number of active elec- 
trons in the n-side generated per unit time is fixed at a given solar 
radiation intensity and cell temperature [9]. Therefore, the current 
curve exhibits a saturation property after the voltage is lower than a 
critical value as shown in Figure 3b. This saturation or choke behav- 
iour is represented properly by the exponential function in Eq. (1). 


2.2 Head-flow rate model for the pumps 


The head-flow rate curves in Figure 2 are very similar to the cur- 
rent-voltage curve of the PV module in Figure 3b. Thus, an expo- 
nential function term should be involved in the mathematical model 
for head-flow rate curves of a partial emission pump. Inspired by the 
mathematical model of PV cells/modules, the following mathemati- 
cal formula is proposed to fit the experimental data in Figure 2: 


H =H, +aQ + bQ2 — ce?2” (2) 


where H and Q are pump head and flow rate, the coefficients 
H,(>0), a, b, c(>0), d(>0), and m(>0) are determined based on 
a set of experimental H and Q data by using the least-squares 
method in Excel, and the number of scattered points of experimental 
data should be more than six. Hg is the head at shut-off point. If 
c=0, then the 2"-order polynomial for fitting head-flow rate curves 
of ordinary centrifugal pumps is recovered. The 2"‘-order polyno- 
mial in Eq. (2) plays a role in fitting the flat part of a head-flow rate, 
but the exponential term takes the dramatic drop in the head-flow 
rate curve into account. 


2.3 Method 


The regression of Eq. (2) was performed against experimental 
data available by employing Solver in Excel. Solver is an add-in 
toolbox. It can be added into Data menu in Excel by clicking File- 
>Options->Add-ins->Solver. The following steps are recom- 
mended to perform a sightly complicated curve fitting with a user 
mathematical model in Excel: 

1) Experimental data of head and flow rate of a partial emission 

pump are read into Excel by occupying two columns; 

2) Assign three columns to accommodate the initial and fitted 
coefficients, statistics results, head predicted with Eq. (2) and 
squared errors between the experimental head and the head 
predicted, respectively; 

3) Calculate the sum of the squared errors and define it as the 
objective function by launching Formula->Sum; 

4) Launch Solver to display the interface for the least-squares 
method by clicking Data->Solver, define the cell for reading 
and showing the objective function, optimization option, op- 
timization algorithm; 

5) Run Solver, if a converged solution has been resulted, check 
the coefficient of determination, R*, compare the experi- 
mental and predicted heads; if not, alter the initial values of 
the fitted coefficients until a satisfactory solution is achieved. 


An example of setup for experimental head-flow rate data, pre- 
dicted head, coefficients, statistical results, and the interface of 
Solver for coefficient optimization is illustrated in Figure 4. In the 
figure, cells H73-H78 are used to store the initial and optimized five 
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Solver Parameters 
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Figure 4. Experimental head-flow rate curve fitting set-up by using Solver in Excel, the right-experimental data, coefficients, model prediction, 
objective function, and statistical results, the left-interface of Solver for coefficient optimization, SS-sum of squares; in Excel cells, SS Total: 


H81=COUNT(D74:D91)*VARP(D74:D91), SS Residual: 


H82=SUMXMY2(D74:D91,E74:E91), 


SS Regression: H83=H81-H82, R?: 


H84=H83/H81, In Excel, COUNT is the function counts the number of cells that contain numbers; VARP is the function that calculates variance 
based on the entire population, SUMXMY7? is the function that returns the sum of squares of differences of corresponding values in two arrays, R? 
is the coefficient of determination, a statistical measure that uses the variance of one variable to explain the variance of another, R°=(SS Total-SS 


Residual)/SS Total 
coefficients, cell H80 is for the sum of errors squared in head, cells 
H81-H84 accommodate the statistical results, cells E74-E91 are uti- 
lized to predict a head with six coefficients and a known flow rate, 
say cell E74 should be assigned with a formula such as 
$H$73+$H$74*C74+ 1 *$H$75*C7442- 
$H$76*EXP($H$77*C74“$H$78) and so on, where cell C74 repre- 
sents the flow rate. If the 2"4-oder term is expected to exclude from 
Eq. (2), the ‘1’ in the formula is altered to ‘0’, ie. 1*$H$75*C74%2- 
>0*$H$75*C74%2. 

In the curve fitting process for Solver, a set of initial values of the 
coefficients Hj, a, b, c, d,and m is needed. A few sets of ini- 
tial values should be tried until Solver is convergent with a better 
coefficient of determination. Once the Solver is convergent with a 
satisfactory coefficient of determination, the coefficients Hp, a, b, 
c, d,and m will be unique. 


3. Results 


Because the experimental data relating to partial emission pumps 
were very limited in the literature, firstly, the head-flow rate curve 
fitting model and method above were applied to five experimental 
partial emission pumps with very low specific speeds suchas n,=23 
[1], n,=20 [10], n,=19 [11], n,=17 [12], and n,=31 [13]. The 
predicted heads are compared with the measurements in Figure 5. 
The determined corresponding coefficients in Eq. (2) are listed in 
Table 1 and 2. The specific speed of partial emission pumps is de- 
fined at best efficiency point or duty point by: 


3.65n,/Q 
n, = 2S5n/e @) 
where n is pump rotative speed, r/min, H is pump head, m, Q is 
pump flow rate, m?/s. The experimental data in Figure Sa, 5d and 5e 
demonstrate the heat-flow rate curves as the impeller rotational 
speed, number of blades and impeller diameter vary. 
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Figure 5. Experimental and fitted head-flow rate curves of existing four par- 
tial emission pumps found in the literature, (a) n,=23 [1], (b) n,=20 [10], 
(c) n,=19 [11], (d) ns=17 [12], (e) n,=31 [13], the lines for fitted curves 
with Eq. (2), Z-number of blades, D2-impeller diameter, the symbols rep- 


resent the experimental data in [1], [10]-[13], respectively 
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Table 1 Coefficients in Eq. (2) for the cases shown in Figure 5a-5c 


n(r/min) 
Coef- : : : 
ficient Figure 5a Figure 5b Figure 5c 
15000 20000 25000 30000 6000 10400 
Ho(m) 178.2723 267.2406 413.2755 567.5908 = 251.4171 421.0318 
a -3.5770 9.2860 1.8542e-13  1.8542e-13 -3.1045 1.1534e-1 
b -29.7348 -1.5236  1.9029e-2.  1.923le-2 9.0563e-1 = -3.0905e-3 
11.2971 1.5490e-4 3.1838e-5 = 8.2953e-5. 1.731 1le-2. 4.3175e-1 
d 5.9866e-5 2.4978e-1 1.0476e-1 1.3744e-1 6.0485e-2 6.0946e-1 
m 46.0748 8.0596 6.8938 5.1233 6.7033 1.0000 
R? 0.8380 0.8882 0.9091 0.8978 0.9954 0.9999 
Table 2 Coefficients in Eq. (2) for the cases in Figure 5d and 5e 
Coeffi- Z in Figure 5d D2(mm) in Figure 5e 
cient 2 4 6 230-216 =~ 202 188 
Ho(m) 76.0433 79.6560 68.0717 84.4175 78.4502 64.9253 58.5929 
a 0.0593 0.3438 0.0682 0.2026 0.8971 0.2248 0.1569 
b ~~ -0.0029 -0.0201 -0.0123 0 0 0 0 
c 0.0418 0.0939 0.0246 0.1888 2.7750 0.0232 0.0174 
d 0.3062 0.6112 0.1736 0.3128 0.3269 1.4085 0.9521 
m 1.0734 0.7420 1.3341 0.8666 0.6648 0.5094 0.6960 
R? 0.9995 0.9926 0.9990 0.9934 0.9930 0.9171 0.9857 


In Figure 5a, the curve fitting is sightly poor because the coeffi- 
cients of determination are R?=0.84, 0.89, 0.91, 0.90 at 15000, 20000, 
25000, 30000r/min, respectively, in Table 1. R? is the coefficient of 
determination reflecting the proportion of the variation in head that 
is predictable from flow rate. In the other cases, as shown in Figure 
5b-e, the curve fittings are quite good since the coefficients of deter- 
mination are as great as R7=0.92-0.99 shown in Table 2. These facts 
suggest that the mathematical model expressed by Eq. (2) is appli- 
cable to the cases with fixed impeller rotational speed and geometry 
but also is adaptive to the change in either impeller speed or geome- 
try. 

Secondly, the experimental data of additional two partial emission 
pumps with a bit high specific speeds at n,=69, 80 in [2] were fitted 
by using Eq. (2) and are compared with the fitting curves in Figure 
6, meanwhile the corresponding fitting coefficients are listed in Ta- 
ble 3. The variation of the experimental head against the flow rate is 
more complicated than that in Figure 5, the fitting is not perfect near 
the region where the head starts to drop off sharply. In other words, 
the fitting curves near that region are smoother than those in the ex- 
perimental data. Consequently, the R? values are in the range of 0.94- 
0.98. Overall, the fitting results are satisfactory. 


4. Discussion 


As shown in Figure 5, the head curve of three partial emission 
pumps is flat in a wide range of flow rate and then drops off with 
increasing flow rate. One may think that the flat part in the head 
curve would be the most suitable for the operation of partial emission 
pumps and it is sufficient to fit the flat part in the curve rather than 
the whole curve. Unfortunately, the pump efficiency in the flat part 


of the curve is lower than in the elbow of the curve, see Figure 7. To 
achieve a better efficiency, a partial emission pump can operate in 
the elbow of the curve. In this case, fitting the whole head curve of 
the pump is necessary. 

The rapid drop in the head with increasing flow rate shown in Fig- 
ure 2 is due to the saturation or choke effect of flow through the dif- 
fuser throat of the volute at a high flow rate. Also, this effect results 
in a large amount of hydraulic loss in the volute or even cavitation at 
the volute tongue [13] and leads to a remarked head drop. Obviously, 
this effect resembles the current saturation or choke behaviour in a 
PV cell. The exponential function of the last term in Eq. (2) is em- 
ployed to handle this effect mathematically. 
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Figure 6. Experimental and fitted head-flow rate curves of existing two 
partial emission pumps found in [2], (a) n,=69, (b) n,=80, the lines for 
fitted curves with Eq. (2), the symbols represent the experimental data in 
[2] 
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Table 3 Coefficients in Eq. (2) for the cases shown in Figure 6 
n(r/min) 
Coefficient Figure 6a Figure 6b 
11000 13000 17000 19000 7200 9500 12000 
Ho(m) 60.5195 92.9577 137.9416 175.0071 42.4756 67.0659 109.5380 

a 3.1587 3.4674 9.2794 7.4943 1.3438 4.7721 2.1250 

b -2.1351 -1.8603 -2.4424 -1.8765 -8.8415e-1 — -1.2698 -5.1330e-1 

c 1.9932e-5 1.5943e-3.  4.5343e-4 2.1998e-5 1.6318e-5 = 2.9735e-4 —-1.4335e-3 

d 4.2687e-3 7.3030e-4 1.0976e-3 3.2434e-3 4.5655e-2 1425le-3  2.1222e-3 
m 6.2347 6.2700 5.2426 4.5320 3.2659 4.9244 4.1350 
R? 0.9459 0.9536 0.9467 0.9623 0.9779 0.9768 0.9866 


260 


6000r/min wit inducer 


230 


Best efficiency point 


0 0.5 2 2.5 


1 15 
Q(m*/h) 
Figure 7. Head H, shatt-power P and etticiency n curves of the partial 
emission pumps tested in [10]. 


In pump industry and its application sectors, the head-flow rate 
curve of centrifugal pumps is usually expressed by using a 2™ or 
even higher order polynomial, for example, in the work presented in 
[14-18]. Furthermore, orthogonal polynomials are employed to fit 
experimental data of centrifugal or axial-flow pumps [19]. Interest- 
ingly, the following mathematical model was proposed to fit the ex- 
perimental data of head of centrifugal pumps [20]: 

H = Hy, —ce®? (4) 
where Hy, c and d are fitting coefficients. Obviously, Eq. (4) is a 
simplified version of Eq. (2) when the conditions: a=b=0 and m=1 
are held. Therefore, the work presented in the paper is innovative to 
partial emission pumps. 

In Figure 5c-5e, the head reduces more gently with increasing 
flow rate than in Figure 5a and 5b, and the head-flow rates resemble 
to the curves of centrifugal pumps, thus the m values are close to 1, 
as shown in Table | and 2. Hence, the mathematical model Eq. (2) are 
suitable for centrifugal pumps as well. 


5. Conclusions 


A mathematical model was proposed for the experimental head- 
flow rate curve fitting of partial emission pumps and implemented 
in Excel. The physics behind the model was explained based on the 
flow rate saturation effect in the pumps. The model can produce 
high-quality curve fitting based on the experimental head-flow rate 
data of seven partial emission pumps collected in the literature and 
is more general than the existing mathematical models. It is hopeful 


that the model can find applications to experimental data processing 
in the pump industry soon. The author can provide the supplemen- 
tary data in Excel file on request. 
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Nomenclature 


a coefficient in Eq. (2) 
b coefficient in Eq. (2) 
Cc coefficient in Eq. (2) 
d coefficient in Eq. (2) 


dD, inlet diameter of impeller (mm) 

D, outlet diameter of impeller (mm) 

H pump head (m) 

Ho coefficient in eq. (2) (m) 

I output current of flat photovoltaic cell/module 
(A) 

Ty diode reversal saturation current shown in Figure 
3a (A) 

lyn photocurrent shown in Figure 3a (A) 

k Boltzmann constant, k=1.38065031*1073J/K 

n rotational speed of impeller (r/min) 

m power in Eq. (2) 

M diode quality factor 

Ns pump specific speed defined by Eq. (3) 

q electron charge, g=1.60217646x10-!°C 

P pump shaft-power (kW) 

Q pump flow rate, (m%/h) or (m*/s) 

R, combined series resistance shown in Figure 3a 
(Q) 

Rsp shunt resistance shown in Figure 3a (Q) 

R? coefficient of determination defined in caption of 
Figure 4 

S solar radiation intensity (W/m?) 

T cell temperature (K) 

V Output voltage of flat photovoltaic cell/module 
(V) 

Z number of blades 

n pump efficiency (%) 
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w rotational angular speed of impeller (rad/s) 
PV photovoltaic 
SS sum of squares 
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